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We report an experimental implementation of a single-qubit generalised measurement scenario 
(POVM) based on a quantum walk model. The qubit is encoded in a single-photon polarisation. The 
photon performs a quantum walk on an array of optical elements, where the polarisation-dependent 
translation is performed via birefringent beam displacers and a change of the polarisation is imple¬ 
mented with the help of wave-plates. We implement: (i) Trine-POVM, i.e., the POVM elements 
uniformly distributed on an equatorial plane of the Bloch sphere; (ii) Symmetric-Informationally- 
Complete (SIC) POVM; and (hi) Unambiguous Discrimination of two non-orthogonal qubit states. 

PACS numbers: 03.65.Ta, 03.67.Ac, 03.67.Lx, 42.25.Hz 


I. INTRODUCTION 

The basic unit of quantum information is a two-level 
quantum system commonly known as a qubit. Qubits 
can be implemented on physical objects such as polari¬ 
sation of photons or intrinsic angular momentum (spin 
1/2) of quantum particles. Any quantum computation 
relies on precise preparations, transformations and mea¬ 
surements of such systems. Before actual quantum com¬ 
puter is build, one has to master the ability to manipulate 
with single qubits and learn how to readout information 
encoded in them. 

The information readout from a quantum system is 
done via a measurement. In most common scenario one 
performs the von Neumann measurement that projects a 
state of the qubit onto one of two perfectly distinguish¬ 
able (orthogonal) physical states of the system. Such 
measurements are sharp in a sense that once the mea¬ 
surement is done, the outcome of the measurement is 
determined and any repetition of exactly the same mea¬ 
surement would yield the same outcome value. 

Physically, von Neumann measurements are realised 
via interaction of the system with the measurement appa¬ 
ratus. The pointer of the measurement apparatus is rep¬ 
resented via wave packed and the interaction causes this 
wave packet to move either to the left or right, depending 
on the value of the measured observable. In general, this 
value might be undetermined and the pointer goes into 
superposition of being to the left and to the right from 
its initial position. The actual collapse of wave function 
is usually attributed to the observer who reads out the 
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measurement outcome from the pointer. The sharpness 
of the measurement comes from the fact that the initial 
spread of the pointer’s state is assumed to be narrow 
and the translation caused by the interaction with the 
measured observable is large enough to prevent overlap 
between the part of the wave packed that was shifted to 
the right with the part that was shifted to the left. 

On the other hand, there are generalised measurement 
scenarios, the so called Positive-Operator Valued Mea¬ 
sure (POVM), in which one allows the principal system to 
interact with an ancillary system, whose state is known, 
and later performs a von Neumann measurement on the 
joint system. This effectively extends the dimensionality 
of the Hilbert space and one can implement measure¬ 
ments of a quibit with more than two outcomes. As a 
result, one gains a plethora of new possibilities. They al¬ 
low one, for example, to perform a tomography of qubit 
with a single measurement setup P, [3| j or to discriminate 
between non-orthogonal quantum states Si- POVMs 
were implemented in laboratories using various setups 

iQ. 

In (T^ it was proposed that they can be implemented 
via a discrete-time quantum walk which has been realised 
in a laboratory using various physical systems [l3l-l^. 
Quantum walks model an evolution of a particle in a dis¬ 
crete space. The particle moves either one step to the left 
or right, depending on a state of a two-level system known 
as a coin. Quantum walks were originally proposed as an 
interference process resulting from a modified version of 
a von Neumann measurement in which a pointer state 
distribution is much broader than the shift of its posi¬ 
tion [ 2 ^ . The position of the pointer plays the role of the 
quantum walker and the qubit that is measured plays the 
role of the coin. If the interaction between the pointer 
and the qubit occurs at the same time as the evolution 
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of the qubit, the measured value changes during the pro¬ 
cess and the pointer starts to move back and forth. This 
movement leads to an interference and the interference 
pattern produced in this process can be interpreted as 
a POVM. In fact, in [l^ it was shown that any POVM 
can be implemented in such a way, provided a necessary 
evolution is applied to the qubit. 

In this work we report an experimental implementa¬ 
tion of the above quantum walk POVM scenario. Here, 
we use the optical setup in which the qubit is encoded 
in a polarization state of a single photon and the posi¬ 
tion of the quantum walker is implemented on a pho¬ 
tonic path 1^. We construct setups realising (i) three 
POVM elements symmetrically distributed on an equa¬ 
torial plane of the Bloch sphere (Trine-POVM); (ii) 
Symmetric-Informationally-Complete (SIC) POVM; and 
(iii) Unambiguous Discrimination of two non-orthogonal 
quantum states. 

II. DISCRETE-TIME QUANTUM WALKS 


III. POVM IMPLEMENTATION WITH 
QUANTUM WALKS 

The probability distribution of a quantum walk, that 
is initially localised at the origin cc = 0, depends on the 
initial coin state a\ —>■) -I- /3| •<— ) and on the subsequent 
coin operations. A single application of the operator T 
makes the particle to go into superposition a|l,—>■) -|- 
/3| — l,t—). In this case the position measurement would 
correspond to a von Neumann measurement of the coin 
in the basis {^,—>•}, where the result —>■ corresponds to 
finding the particle at a; = 1 and •<— to finding it at a; = 
— 1. However, if one allows quantum walk to evolve for 
more than one step and the coin operation to change from 
one step to another step, then the particle can spread 
over many positions and the measurement of its location 
at X may correspond to a measurement of some POVM 
element Ex- Indeed, a proper choice of coin o pera tions 
can lead to an arbitrary qubit POVM scenario |1^ . 

The POVM elements Ei (i = 1, 2,..., n) are nonnega¬ 
tive operators obeying 


Discrete-time quantum walks are quantum counter¬ 
parts of classical random walks in which a particle takes 
a random step to the left or right. In the classical case 
the particle spreads in a diffusive manner (a standard 
deviation of its position is proportional to \/t) and after 
many steps a spatial probability distribution is Gaussian. 
In quantum case the system is described by two degrees 
of freedom \il)) = \x, c): the position a; = ..., —1,0,1,... 
and a two-level internal degree of freedom known as a 
coin c =<—,—>■. The evolution is unitary and consists of 
two sub-operations U = CT^ namely a unitary coin toss 
(that is a 2 X 2 unitary matrix acting only on the coin 
degree of freedom) 

/ cos 9 e“*^sin0\ . . 

y—e*^sin6* cos0 j 

and a conditional translation operation 


Y,E, = 1. (4) 

i 

They differ from standard von Neumann projectors H^ in 
that they do not have to be orthonormal (HiH^- = SijTli, 
but EiEj ^ 6ij) and that their number can be greater 
than the dimension of the system n > d. The follow¬ 
ing quantum walk algorithm, proposed in [l^ , generates 
an arbitrary set of rank 1 POVM elements {Ei ,..., En} 
(rank 2 elements can be generated with a modified ver¬ 
sion of this algorithm). 

1. Initiate the quantum walk at position a; = 0 with 
the coin state corresponding to the qubit state one 
wants to measure. 

2. Set i := 1. 


3. While z < n do the following: 


T — |a: -I- 1, —>■ I -|- |x — 1, ■^){x, t— |. (2) 

X 


Quantum walks with uniform (position-independent) 
coin operation spread ballistically (standard deviation 
proportional to t) and its probability distribution differs 
form the classical Gaussian shape (^ . . On the other 

hand, quantum walks with position-dependent coin op¬ 
eration 


Cx 


f cos 9x 



( 3 ) 


• Apply coin operation at position a; = 0 
and identity elsewhere and then apply trans¬ 
lation operator T. 

(o') 

• Apply coin operation C) ^ at position x = 1, 


NOT = 



at position x = —1 and identity elsewhere and 
then apply translation operator T. 

• i ■= i + 1. 


can be used to observe localisation [3l|,[s^ , or to simulate 
physical systems with non-homogenous interactions . 
In [l^ it was shown that quantum walks with position 
and time-dependent coin operations Cx,t can be also used 
to implement POVMs. 


Since in this work we are only interested in the measure¬ 
ment statistics, and not in the pos t-measurement states, 
we simplified the algorithm in and omitted the last 
step. The POVM elements that are generated depend 
solely on the form of operators and 
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IV. EXPERIMENT 


A. Trine POVM 


In our experiment, frequency-doubled femtosecond 
pulse (390 nm, 76 MHz repetition rate, 80 mw aver¬ 
age power) from a mode-locked Ti:sapphire laser pump a 
type-I beta-barium-borate (BBO, 6.0x6.0x2.0 mm^, 6= 
29.9) crystal to produces the degenerate photon pairs. 
After being redirected by the mirrors (Ml and M2, as 
in the Fig. 1(a)) and the interference filters (IF, AA=3 
nm, A=780 nm), the photon pairs generated in the spon¬ 
taneous parametric down-conversion (SPDC) process are 
coupled into single-mode fibers separately. Single photon 
state is prepared by triggering on one of these two pho¬ 
tons, and the coincidence counting rate collected by the 
avalanche photo-diodes (APD) are about 4 x 10“* in one 
minute. 



Coincidence 


Dectection 




The experimental setup in Fig. 1(a) is used to con¬ 
struct the trine POVMs, = 15 2,3), 


1^3) = \ H ) 

\^l) = -lm-V3\V)) (5) 

m = -lm+V3\v)). 


According to the settings of the coin operators, the opti¬ 
cal axes of BDl and BD2 must be aligned, in other words, 
they form an interferometer. When rotating HWPl and 
HWP3 by 22.5°, we observed that the interference visibil¬ 
ity of the interferometer was about 99.8 and the system 
was stable over 2.5 hours of timescale. After aligning, 
we begin to set the corresponding coin operators in each 
step. For the Trine POVMs, we have 




1 , 0 ) 


( 2 ) 






( 6 ) 


FIG. 1: Experimental setup, (a) Experimental setup for con¬ 
structing the trine POVMs corresponding to \ipl) and realiz¬ 
ing the unambiguous state discrimination, (b) Optical net¬ 
work for constructing SIC-POVMs. Initial coin states are 
prepared by passing the single photons through a polarizing 
beam splitter (PBS), a half-wave plate (HWP) HWPl and a 
quarter-wave plate (QWP) QWPl in a specific configuration. 
The conditional position shifts are implemented by Beam Dis¬ 
placers (BDs) and the coin operators in different positions are 
realized by wave plates with different angles (Table 1). The 
indices in the figure denote the position of walker. 


One-dimensional discrete time quantum walk system 
has two degrees of freedom, x and c, x is the position 
of the particle and c is the state of the coin. In our 
experiment they are encoded in the longitudinal spa¬ 
tial modes and polarizations \H), \V) of the single pho¬ 
tons respectively. In this case, the conditional transla¬ 
tion operator as given by Eq.(2) is realized by the de¬ 
signed BD, that does not displace the vertical polarized 
photons {\x,V) —>■ |a; — 1,F)) but makes the horizon¬ 
tal polarized ones undergo a 4 mm lateral displacement 
(\x,H) -)> \x + l,H)). 


where and are realized by rotating HWP2 and 
HWP3 by 17.32° and 22.5° respectively. The initial trine 
coin states |'03) are constructed by rotating HWPl by 
0°, —30° and 30°, while the anti-trine states j'i/j^) are 
constructed by rotating it by 45°, 15° and —15°. At last, 
every output port’s detection efficiency are calibrated so 
that the differences among them are below 5%. 

Fig. 2 shows that the results in our experiment agree 
with theoretical predictions. The ratios 2/3 : 1/6 : 1/6 
(0:1/2: 1/2) for the cases of l^g) (I'l/g), where = 

0) are given in theory and the detailed numerical results 
of the probability distributions can be found in Table H 
and HI. To visualize that the setup has constructed the 
Trine POVMs, it is important to demonstrate that we 
cannot find states I'i/jl) in position 4, iVil) in position 0 
and I'i/I) in position 2. Fig. 2(b) and Table HI show that 
the probabilities of these events are indeed very close to 
zero, with an average value of 0.0085. In addition, the 
results for states \tjjl) also indicate the coefficients of the 
POVM we constructed is |, see Fig. 2(a). The errors 
in our experiment mainly stem from the imperfect wave 
plates and the interferometers and the counting statistics 
of the photons. 
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FIG. 2: Results for trine POVMs. Histogram shows the prob¬ 
abilities of counting rates in position 0, 2 and 4 with input 
states tpKa.) and tpl{h), respectively. All results are normal¬ 
ized so that the sum of the counts in these three positions is 
1. The theoretical values are shown as the blue lines, which 
are 2/3, 1/6, 1/6 for ipl and 1/2, 1/2, 0 for 2, 3); 

error bars are too small to identify. 


B. SIC POVM 


The optical network in Fig. 1(b) construct the SIC 
POVMs, i|V^4)(V-l|(t = l,2,3,4), 


IV'i) = \H) 

\^l) = -±\H) + e^^^\V) 
The coin operators 


cf ^ = 

C'« = 


1,C 

1 

V2 

1 

71 


( 2 ) ^ J _ (-1 



73U 

,cP = i 



( 8 ) 


are realized by wave plates in various configurations (de¬ 
tails in Table I). 



ct 

ck 

ci 


ci 

HWP 

67° 30' 

67° 30' 

17° 38' 

142°30' 

— 

QWP 




150° 

— 


TABLE I: The configurations of the QWPs and HWPs to 
realize the coin operators for constructing the SIC-POVM. 


For these settings, BDl and BD2, BD3 and BD4 form 
two interferometers whose interference visibility are both 
above 0.993. The HWPl and QWPl with different angles 
in front of a polarizing beam splitter (PBS) are used to 
produce '04 and the corresponding orthogonal states 'tjj\ 
(Table VI and Table VII). As shown in Fig. 3, the results 
are also in good accordance with theoretical ratios ^ ^ ■ 

i : i for 04 and ^ ^ ^ : 0 for 0 |. 


(a) 


(b) 





^^IpositiorO 
I I positior? 
1^ < positior4 
^^BpositiorS 


llllllftiBI 






FIG. 3: Results for SIC POVMs. Histogram showing the 
normalized probability of counting rate in position 0, 2, 4 
and 6 respectively with the input state 04(a) and 04(b). The 
theoretical values are shown as the blue lines, which are given 
by 1/2, 1/6, 1/6, 1/6 for states 04 and 1/3, 1/3, 1/3, 0 for 
states 04; error bars are too small to identify. 


C. Unambiguous state discrimination 


For the unambiguous state discrimination of states, 
| 0 ±) = cos{9/2)\H) ± sin(0/2)|V), we can use the same 
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on a quantum walk model presented in (l^ . Our results 
match the theoretical predictions. We believe that these 
kind of experimental setups can be used in the future to 
implement other types of generalised measurement sce¬ 
narios with multiple outcomes and rank 2 POVM ele¬ 
ments, and to study quantum walks with position and 
time-dependent coin operations. Finally, we would like 
to mention that similar results had been reported while 
we were working on this experiment . 


VI. ACKNOWLEDGEMENTS 


FIG. 4: Theoretical and experimental successful probability 
v.s. 9 which is related to the state to be discrimination; error 
bars are too small to identify. 


setup as in Fig. 1(a), with 

Cf) = 1, 


cf ^ = 


^1 — (tan |)2 tan | 

tan I — — (tan 


( 9 ) 


The state after four quantum walk steps becomes 

0 . 


|'0+) = Vcos6\4:,H) + y/2sin-|2, iJ), (10) 


Q 

= \/cosd\4:,H) — y/^sin-|0, iJ). (11) 

Therefor, if the photon is detected at position x = 0 one 
knows that the coin was definitely in the state \ip-). If 
it was detected at position x = 2 one knows that the 
coin was definitely in the state |'0+). Finally, if it was 
detected at position x = 4 one gains no information and 
the discrimination was unsuccessful. 

In our experiment, the input states and ip- are 
prepared in ^ steps from 0 < 0 < . The states for var¬ 

ious 9 are prepared by rotating the HWP, placed before 
the polarizing beam splitter (PBS), and the coin oper¬ 
ator c[^\ C^^^^and the NOT operator are realized by a 
HWP rotating in the angle of i arcsin(tan |), ^ and j 
respectively. From Fig. 4 we can see that the probability 
of the successful discrimination is increasing with 9. For 
more details see Tables VIII and IX. 


V. CONCLUSIONS 
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VII. APPENDIX 


The detailed results about tpl, Pp\, ip\ and 'ipX are shown 
in Table II, Table III, Table IV and Table V. Table VI 
and Table VII are the angles of the HWPl and QWPl for 
preparing states tp\ 'ip\. The main parameters of un- 
ambignous state discrimination and the detailed results 
can be found in Table VIII and Table IX. 


state 

Po 

P2 

Pi 

tpz 

0.1684(20) 

0.1711(19) 

0.6604(25) 

V's 

0.6540(26) 

0.1731(20) 

0.1729(21) 

1*3 

0.1753(21) 

0.6466(25) 

0.1782(21) 


TABLE II: Po, ft and P 4 are the nomarized probabilities of count¬ 
ing rate in position 0, 2 and 4. 


state 

Po 

P2 

Pi 

V’s 

0.5018(27) 

0.4977(27) 

0.0005(01) 

V'l 

0.0151(06) 

0.4897(27) 

0.4952(26) 


0.4933(25) 

0.0081(05) 

0.4987(25) 


We experimentally realised three generalised measure¬ 
ment scenarios for a qubit. These scenarios are based 


TABLE III: Po, P 2 and P 4 are the nomarized probabilities of 
counting rate in position 0, 2 and 4. 
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state 

Po 

P 2 

Pa 

Re 

ipl 

0.1662(19) 

0.1654(19) 

0.1934(20) 

0.4749(26) 

i’l 

0.1585(18) 

0.1571(18) 

0.5220(27) 

0.1625(20) 

V'l 

0.5015(26) 

0.1676(19) 

0.1695(19) 

0.1614(19) 

V'l 

0.1885(20) 

0.4843(25) 

0.1623(19) 

0.1649(19) 


TABLE IV: The nomarized probabilities of counting rates for state 
r/)! in position 0, 2, 4 and 6 . 


state 

Pq 

P 2 

Pa 

Pe 

ipl 

0.3341(24) 

0.3367(25) 

0.3289(24) 

0.0003(01) 

ipl 

0.3182(23) 

0.3283(24) 

0.0051(04) 

0.3485(24) 


0.0040(03) 

0.3209(24) 

0.3671(25) 

0.3080(23) 


0.3152(24) 

0.0005(01) 

0.3647(24) 

0.3196(24) 


TABLE V: The nomarized probabilities of counting rates for state 
i/il in position 0, 2, 4 and 6 . 



ipl 



V'l 

HWPl 

0° 

-27° 22 ' 

17° 38' 

45° 

QWPl 

0° 

35° 16' 

-27°22' 

-27° 22 ' 


TABLE VI: The angles of HWPl and QWPl used to prepare 
the states 



ipl 


V>4 


HWPl 

45° 

17°38' 

-27° 22 ' 

0° 

QWPl 

0° 

35°16' 

-27° 22 ' 

- 2 T 22 


TABLE VII: The angles of HWPl and QWPl used to prepare 
the states 


e 

Sl/2 


Ptheory 

Pe 

n /20 

2° 15 ' 

2° 15 ' 

0.0123 

0.013(06) 

tt/IO 

4°30' 

4° 34 ' 

0.0489 

0.050(10) 

37 r/20 

6° 45 ' 

6° 57 ' 

0.0109 

0.103(15) 

7 r/5 

9° 

9° 29 ' 

0.191 

0.181(19) 

7 r/4 

11°15' 

12° 14 ' 

0.293 

0.285(22) 

37 r/10 

13°30' 

15° 19 ' 

0.412 

0.402(23) 

77 r/20 

15°45' 

18° 54 ' 

0.546 

0.531(24) 

2n/5 

18° 

23°18' 

0.691 

0.673(22) 

97 r/20 

20°15' 

29°2l' 

0.844 

0.832(18) 

7r/2 

22°30' 

45° 

1.000 

0.996(03) 


TABLE VIII: 9, the angle related to the input state which 
is prepared by HWPl rotated to 61 / 2 ; Cf, the angle of HWP2 to 
realize the operator Cj. Pthnory and Pg represent the theoretical 
and the experimental successful probability. 
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Sl/2 


Ptheory 

Pe 

-7r/20 

-2°15' 

2° 15 ' 

0.0123 

0.0127(05) 

-tt/IO 

-4°30' 

4° 34 ' 

0.0489 

0.0511(11) 

-37r/20 

-6°45' 

6° 57 ' 

0.0109 

0.1096(15) 

—7r/5 

-9° 

9° 29 ' 

0.191 

0.1889(19) 

—7r/4 

-11°15' 

12° 14 ' 

0.293 

0.2883(23) 

-37r/10 

-13°30' 

15° 19 ' 

0.412 

0.4092(24) 

-77r/20 

-15°45' 

18° 54 ' 

0.546 

0.5340(24) 

-27r/5 

-18° 

23°18' 

0.691 

0.6812(22) 

-97r/20 

-20° 15 ' 

29°2l' 

0.844 

0.8367(17) 

-7r/2 

-22°30' 

45° 

1.000 

0.9951(03) 


TABLE IX: 9, the angle related to the input states ip— which is 
prepared by HWPl rotated to 9i/2', Cf, the angle of HWP2 to 
realise the operator Cf. Ptheory and Pe represent the theoretical 
and the experimental successful probability. 
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